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We study sequential tunneling of magnetic excitations in nonitinerant systems (either magnons 
or spinons) through triangular molecular magnets. It is known that the quantum state of such 
molecular magnets can be controlled by application of an electric- or a magnetic field. Here, we use 
this fact to control the flow of a spin current through the molecular magnet by electric- or magnetic 
means. This allows us to design a system that behaves as a magnon-transistor. We show how to 
combine three magnon-transistors to form a NAND-gate, and give several possible realizations of 
the latter, one of which could function at room temperature using transistors with a 11 ns switching 
time. 

PACS numbers: 75.50.Xx, 75.30.Ds, 75.76.+j 



I. INTRODUCTION 

In spintronic devices in insulating magnets, informa- 
tion about the logic state can be encoded in collec- 
tive magnetic excitations, typically either spinons or 
magnons.— Due to the nature of these carriers, power 
dissipation in such nonitinerant devices is anticipated to 
be much lower^ than in their electronic counterparts, 
as well as in spintronics devices in semiconductors^ As 
excess heating is a limiting factor in modern electron- 
ics, spintronics in insulating magnets is considered a con- 
tender to become the next computing paradigm^— 

Since any classical algorithm can be implemented us- 
ing a combination of transistors only, the design of this 
element in insulating magnets is a pivotal issue. Here, we 
theoretically show that that it is possible to make a tran- 
sistor in which the logic state is encoded in purely mag- 
netic excitations, and whose operation can be controlled 
by either a magnetic- or electric field. In our transis- 
tor, triangular molecular magnets^— take the role of 
gate, and we model the source and drain by spin reser- 
voirs. We show that our transistor, which could operate 
at high clock speed at room temperature, can be used 
to implement the NAND-gate, one of the two existing 
universal gates for classical computation. 

In our system, transport of magnetization occurs pri- 
marily by sequential tunneling of magnons (for FM reser- 
voirs) or spinons (for AF reservoirs) through the molec- 
ular magnet. We will show how it is possible to sup- 
press or increase this sequential tunneling (and thereby 
switch between the insulating and conducting state of 
the transistor) by changing the internal state of the 
molecule through external fields, either electric or mag- 
netic. Molecular magnets are especially suitable due to 
their chemical variety and controllability, as well as their 
relatively large size, which makes control of the state 
easier. For similar reasons, they have been proposed as 
good building blocks for novel spin-polarized—, as well 
as quantum computing devices i^r— 

This work is organized as follows: In Sec. [IT] we intro- 
duce in more detail the previously mentioned system in 
which we will create our transistor. In Sec. IIIII we cal- 



culate the tunneling rates of magnons / spinons through a 
triangular molecular magnet, and calculate the spin cur- 
rent through the molecular magnet from these rates. In 
Sec. II VI we show how controlling the state of a molecular 
magnet by electric- or magnetic fields allows us to design 
a transistor for either magnons or spinons. In Sec. [V] 
we focus on possible implementations of our transistor. 
Finally, we discuss certain constraints on our results in 
Sec. ED 



II. SYSTEM 

The system we employ consists of a triangular molec- 
ular magnet, such as {CU3} (see Ref. |2l[ ), with initially 
two vertices weakly exchange-coupled to identical spin 
reservoirs, see Fig. []Ia)-(b). We will consider both 
one-dimensional (ID) and two-dimensional (2D) FM spin 
reservoirs as well as ID AF spin reservoirs. 

Both spin reservoirs, the molecular magnet, and the 
weak coupling between the subsystems are initially de- 
scribed by the isotropic Heisenberg Hamiltonian with 
Dzyaloshinkii-Moriya (DM) interaction 

H = J2 -'-A ' S j + D • {St x Sj) . (1) 

(ij) 

The exchange interaction Jy and DM vector = Dije z 
are assumed to be constant for each subsystem, with val- 
ues given in Fig. QJb). We will assume that the J"s set 
the smallest energy scale in the system, so that we can 
analyze tunneling processes using perturbation theory. 
In our model, a finite spin current is induced by appli- 
cation of a magnetic field AB to the left spin reservoir, 
which creates a non-equilibrium distribution of magnetic 
excitations. In reality, due to the finite lifetime of the 
magnetic excitations, a steady state spin current has to 
be generated using e.g. an AC magnet field difference, a 
static temperature difference, or spin pumping. 

When Dm *C Jm [see Fig. QJb)], the low-energy sub- 
space of the triangular magnetic molecule consists of a 
quadruplet with total spin- 1/2, and the eigenstates of the 
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FIG. 1. (a)-(b) Pictorial representation of a single purely magnetic spin transistor, including Heisenberg parameters of the 
subsystems and the different magnetic- and electric fields. The field Bq is applied to both spin reservoirs (here shown as f D 
spin chains) and the molecular magnet, (c)-(d) Excitation spectrum corresponding to Eq. ((2j), which can be rewritten in 
diagonal form as Hm = gml-iBBSg + v/ (Dm /2) 2 + (dE') 2 C'g, where the subscript denotes an operator in the basis given by 
| t> l)fl — sin 6*| t, f ) + cos 6\ f, — 1), | t> — 1) = — cos 6\ t, 1) + sin6!| — 1), and | |, | 4-j — )e the same but with the spins 
flipped. We defined tanfl = [^(D M /2) 2 + {dE') 2 + D M /2]/dE'. We have absorbed the spin-orbit interaction in the definition 
of the basis by performing a spin-dependent rotation in chirality space. In (c) we choose gM/xs-B = Dm — 10 -3 , in (d) we put 
L>M = 5-I(r 3 . 



Hamiltonian can be labeled as \ms,nic}- These states 
are eigenstates of the z projection of two mutually com- 
muting effective spin-1/2 degrees of freedom: the total 
spin operator S (with eigenvalues mg = ±1/2) and the 
chirality operator C (with eigenvalues mc = ±1). An in 
plane electric field E = E x e x + E y e y couples to the chi- 
rality through the spin-electric effect ; 22 ' 23 the coupling of 
a magnetic field B = Be z to the total spin is described 
by the Zeeman interaction. The Hamiltonian for the low- 
energy subspace is 

H M = 9M[i B BS z + dE' ■ C" + D M S Z C Z . (2) 

Here, B = B' + Bq, d is the effective dipole moment 
of the molecule, E 1 is the rotated electric field^ 2 - and 
C" = C x e x +C y e y . The g- factor of the molecular magnet 
is denoted gu- Diagonalization of Eq. ^ for E' = E'e x 
gives the spectrum depicted in Fig. [Hc)-(d). 

Tunneling of magnetic excitations between the left- 
/right spin reservoir and the molecular magnet is de- 
scribed by the Heisenberg exchange interaction between 
Sl and Si respectively Sr and S2. By evaluating the 
relevant matrix elements of the Hamiltonian in the low- 
energy subspace of the molecular magnet, we find the 
respective effective Hamiltonians 

Hlm = y [1 - 2sin(20)C| + 2 cos(20)C£ ] S L • S 9 , 



Umr — —r 
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1 + sin(26»)Q - cos(260C£ + V3C l e 



Se • Sr. 

(3) 

We interpret Eqs. ([3]) as effectively describing tunneling 
of spin excitations from the spin reservoir onto the total 



spin of the molecular magnet and vice versa, with a tun- 
neling strength that depends on the chirality state of the 
molecule. 

We will also consider the case in which the third vertex 
of the molecule, described by the spin S3, is coupled to 
a single spin Sc in a separate reservoir. Within the low- 
energy subspace of the molecule, we find the Hamiltonian 



H, 



CM 



1 + sin(20)Cf - cos(20)C£ - VsCj 



■Sc- 

Furthermore, in the case where a single spin S^/r is 
exchange-coupled to all three spins of the molecular mag- 
net with equal strength J', the effective Hamiltonian is 
simply J'S l /r ■ Sg. In this case, the tunneling of mag- 
netic excitations no longer depends on the chirality of the 
molecule. 



The first thing we see from Eqs. ©-(Ill) is that a static 
equilibrium magnetization Se z of a reservoir acts as a 
constant perturbation on the state of the molecular mag- 
net through the relevant exchange interaction between 
reservoir and molecule. We will show that it is possi- 
ble to make this effect trivial, or even beneficial to our 
purposes, in all cases under consideration. Additional 
dynamics of the systems is due to the behavior of mag- 
netic excitations which exist on top of the equilibrium 
magnetization. We will study this dynamics next. 
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III. TRANSITION RATES AND SPIN 
CURRENT 



To determine the spin current through the molecular 
magnet, we use a master-equation approach. We assume 
that energy is conserved in all tunneling processes, and 
ignore higher order effects. Transition rates from ini- 
tial state \i) to final state |/) of the system due to the 
tunneling processes described by Eqs. ([3]) are described 
respectively by Rfj and R^. Using the golden rule, we 



rates 



can calculate Rff R to second order in J[ 2 



as 



R L 1 



dr{i\H' LM (T)\f){f\H' LM (0)\i), (5) 



The same expression holds for Rf* with H' LM replaced by 
H' MR . The apostrophe denotes an operator in the inter- 
action representation. The nontrivial part of the prob- 
lem reduces then to finding correlation functions such as 
(S^, r (t)S^^ r (0)) in the spin reservoirs. In the next two 
sections we will find the relevant expressions for both FM 
and AF spin reservoirs. 

Calculation of the spin current requires both the tran- 
sition rates as well as the probabilities Pi that the 
molecule is in the state \i). We define the vector P = 
(Pf+, Pf-, Pi+, Pi-)- The time evolution of the occupa- 
tion probability vector P is then given by dP/di = RR, 
where R is the 4x4 matrix that contains the appro- 
priate transition rates. The steady state probabilities 
are contained in the kernel of R, normalized such that 
J2i Pi = 1- Hence, P is uniquely determined by the tran- 
sition rates. The spin current Is is then defined as the 
net rate with which excitations leave the left reservoir 
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Next, we will calculate the relevant transition rates for 
the different types of reservoirs. 



FM reservoirs 



R 



T+4-+ 



(7) 
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^fm(wb + uje — was) 



The energy scales are given by Hujb = 9m^bB, Fujab — 
g R jj, B &B, and hu E = 2 v / (D M /2) 2 + (dE') 2 . Also, 
= 1 - 2sin(20) and 77^ = 2cos(20). Furthermore, 
K FM (u) = Pid(^ 9 ) [1 + n B(^ q )}\ aJq=0J , where ft D (w g ) = 

a\duj q /dq\~ is proportional to the density of states 
(DOS) in the spin chain, and ns(uJq) is the Bose-Einstein 
distribution of the magnons. To get the other rates, we 
use that the only effect of inverting the total spin (in 
the initial and final state simultaneously) is to replace 
1 + riB{<-<Jq) — > rising) and oj^ -> —we, and inverting 
both the chiralities changes uje ^ -w^, £+—>•£_ = 
1 + 2sin(20), and 77+ -> r\t = r}\. To obtain the rates 
with respect to the right spin chain, we put was = in 

l±sin(20), 



> ,r 2 . 



Due to en- 



Eqs. 0. Furthermore, we replace ^4 
t]^ -> 77^ = |cos(20) + i-\/3|, and J[ 
ergy conservation, processes that do not flip the total 
spin can only occur for uje = 0. Since we will always 
consider finite lo e , we can put these rates to zero. 

To obtain the rates for a system with a 2D FM reser- 
voir, we simply replace pm(wq) by the 2D DOS, which 
for small |q| is given by /92D(w q ) = h/(AS\J\). 



B. AF reservoirs 

Next, we will derive the transition rates for hopping 
of spinons between semi-infinite AF spin-1/2 chains and 
a triangular molecular magnet. In order to do so, we 
start by giving a description of the spin chains in terms 
of Luttinger liquid theory, which turns out to be a con- 
venient framework for our purpose. For concreteness, 
we focus on the description of the left spin chain. Eq. 
(fTj) can be mapped on the anisotropic Heisenberg Hamil- 
tonian with anisotropy A = J/V J' 2 + D 2 by perform- 
ing a position-dependent rotation in spin space. After 
performing a Jordan- Wigner transformation, taking the 
continuum limit of the resulting fermionic Hamiltonian, 
and subsequent bosonization, the resulting Hamiltonian 
describing the left spin chain [for which x € (—00, 0)] 
becomes 



We first consider the simplest case of ID FM reservoirs 
in the absence of DM interaction, i.e. with D = 0. Us- 
ing the Holstein-Primakoff transformation, we can map 
a ID FM system with S > 1/2 on a system of non- 
interacting bosonic particles (magnons) with dispersion 
huj q = 4| J|S f sin 2 ((7a/2) + gnPBP>o, where gn is the g- 
factor of the reservoir. This allows us to find the required 
correlation functions by rewriting r (t)S^^ R (0)) etc. 

in terms of the bosonic operators a q ,a q . We find the 
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Ax 



uK (d x d{x)f + ± (d^x)) 2 



(8) 



The bosonic density field <fi(x) and its conjugate momen- 
tum field 0(x) satisfy [(f>(x), d x i9{x')} — in5(x — x'). The 
sound velocity u of the bosonic excitations as well as the 
interaction parameter K can be determined from the pa- 
rameters J and D of the spin chain using Bethe Ansatz 
results. At the isotropic point K = 1/2, and K = 1 
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FIG. 2. (a)-(b) log 10 (7s) versus AB for different values of the magnetic field B, for ID AF reservoirs with different Luttinger 
liquid parameters K. Due to the gapless nature of the spinons, we can always set B = Bo and B' = when considering AF 
reservoirs, (c) Idem for ID FM reservoirs. Here, Bo — 1 and B' « B. (d)-(e) log 10 (7s) versus AB for different values 
of dE' , for AF reservoirs with different Luttinger liquid parameters K. We choose Bulks = 0.6 K and B — 150 mT. If we 
assume that the effective dipole moment lies between d — (10 -4 — l)eR (see Ref. [23T I. where R w 1 nm is the bond length of the 
molecular magnet, then dE' /ks = 0.1 K corresponds to an electric field E' ~ (10 8 — 10 4 ) V m . (f) Idem for the FM system. 
Here, Bo = 1 /xT and B' = 150 mT. (g) Illustration of the alternative switching mechanism for FM reservoirs. When the level 
splitting of the molecular magnet is smaller than the minimum energy of a magnon in the lead, the system is in the insulating 
phase. Again, Bo = 1 /iT. The plots (a)-(g) are for parameters J '/ks = 100 K, T — 10 mK, and Ji/fcs = J^/fcs = 0.1 K (see 
text). For the FM plots, 3 = 1. (h) Is versus AB for different values of dE' for the AF system, at an experimentally accessible 
temperature. Parameters are J/k B = 100 K, J' /k B = 2 K, D M /k B = 0.3 K, B = 75 mT, and T — IK. 



describes the free model. To analyze the hopping be- 
tween spin chain and molecular magnet, we will need the 
continuum form of the spin operators in the spin chain. 
These are given by 

i8(x) 

S-(x) = -=(~l) x [cos (2<j>{x) - 2k F x) + 1] , 

V2tto ( 9 ) 

S z (x) = --d x Mx) + — cos \2<f>(x) - 2k F x] . 
7T 7ra 

Here, a is the lattice spacing of the spin chain, and 
kp = n/2a + gupsBo/uh is the Fermi wave vector. For 
simplicity, we assume that Bq = in our derivation. 
Since the spin chain is semi-infinite (with its last spin 
at x = —a), we require that the wave function vanishes 
at the origin. This constrains the density field to a con- 
stant value at the origin, such that cos 0(0) = 0. To 
analyze the behavior of the spin fields at the origin, we 
introduce chiral fields 4>l/r( x ) — KQ(x) ± 4>{x), which 
are related on the entire space by the constraint on the 
density field at the origin^ This allows us to map Eq. 
([8]) on a quadratic Hamiltonian that only depends on 
4>l(x). Performing a renormalization group (RG) analy- 
sis on the spin operators near the boundary then yields 
that S z (x) is a marginal operator, and 5 ± (0) scales as 
1 — 1/(2K), so that it is relevant for systems with finite 
DM interaction. Since the RG flow is stopped either by 
temperature or by the relevant energy scale Em of the 
molecular magnet, this gives the constraint on the hop- 



ping J' [ma,x(k B T, E M ) I J] 



-l + l/(2if) 



•C J for our sequen- 



ces = and Jm,Dm ^> J,D, the entire system can be 
described as two spin chains connected by a single impu- 
rity spin. It is known that, at half-filling (i.e. Bq = 0), 
such a system 'heals' itself at low temperatures^ and 
the current at small AB becomes that of the free sys- 
tem, Is = (g[J,B) 2 AB/h. In this regime, the current in 
the conducting state is strongly enhanced compared to 
that in the sequential tunneling regime. 

In order to calculate the required transition rates, we 
need to calculate spin-spin correlation functions at the 
boundary. Since the density field is frozen out there, the 
sole relevant correlation function is that of the momen- 
tum field 9{t) = 9(0, t). At finite temperature T, it is 
given by 



([9(t) - 9(0)} 



K 



In 



ihijjQ 



sinh 



/tt^o [t-i5] 

\ n 



tial tunneling approach to be valid. We note that when 



(10) 

Here, #o = fesT, and uic is the UV-cutoff of the theory. 
For this model it is approximated as ujq ~ J/h. 6 is 
a positive infinitesimal. The analysis of the right spin 
chains goes along the same lines, and we will refrain from 
repeating the steps here. 

To calculate the transition rates for AF reservoirs we 
substitute Eqs. (0) in Eq. Using the correlation 

function Eq. (fit))) and the fact that <j>(x) is constant at 
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the boundary then gives the rates 

/ j> qL \ 2 
= ^z— - K AF (u B - was), 

R^+l- = [~3j~J k af{u b + we - was)- 

Here, C+ = d L [l - 2sin(20)]/2, v\ = d L cos{29), and 
d,L = \plpn. We have ignored a small fc^-dependent 
contribution to 4 here. The function K\p(ui) describes 
the influence of the spin chain on the transition rate and 
is given by 

/oo 
dre i"T e -i<[0(T)-0(O)] 2 ) 
- OO 

(<*r\ ' 2+ ^ p ™/ UT Wg) + ^/^)| 2 

= " t UJ 6 iW*) 

(12) 

Here, u>t — 2tt8q/H. To get the other rates, we use 
that the only effect of inverting the total spin (in the 
initial and final state simultaneously) is to change the 
sign of lob and was, and inverting both the chiralities 
changes uj e -> -Ub, C+ -> C- = ^[1 + 2 sin(20)]/2, and 
— ^ i/f; = z/f. To obtain the rates with the respect 
to the right spin chain, we put was = in Eqs. (TTT]) . 
Furthermore, we replace C± — » C± = <^r[1 ± sin(20)]/2, 
z/£ -> v£ = d R \cos(29) + iV3\/2, and J[ -> J£. The 
constant dp = d^. For the AF reservoirs, as for the FM 
reservoirs, processes that do not flip the total spin can 
only occur for uje — 0. Since we will always consider 
finite lje, we can put these rates to zero. 

At this point, we have derived the transition rates 
for both AF and FM reservoirs. As we have seen, the 
main difference in the resulting rates is the replacement 
of the bosonic DOS and distribution function by power 
law functions, typical for Luttinger liquid models. We 
note here that, at low energies, the bosonic character of 
the magnons in the FM system yields larger spin currents 
than the fermionic spinons, which is extremely beneficial 
for the application we have in mind here. 



be switched between the insulating- (for ujb ~ wb 2 ) and 
conducting (for lub ~ ^>b 2 ) state. This is shown in Fig. 
Ha)-(c). 

The system with FM reservoirs offer an additional pos- 
sibility to switch between the insulating and conducting 
state: When was > ojb, the minimum energy of the 
magnons in the reservoir exceeds the level splitting of 
the molecular magnet. In this case, the system is also 
insulating (neglecting higher order processes). This be- 
havior has been indicated in Fig. [2jg). The use of this 
mechanism to switch between insulating and conducting 
states requires smaller magnetic fields compared to the 
previously discussed mechanism. 

Alternatively, the switching behavior can be controlled 
by an electric field. The mechanism is different from 
that for magnetic control, since an electric field does not 
increase the splitting between the two lowest states with 
opposite total spin. However, we note that when fcgT <C 
hwE, transport occurs through transitions between states 
in the subspace spanned by | f, — l)g and | J,, — l)g. Eqs. 
([3]) show that when 8 — > 7r/4, i.e. when dE' 3> D m, 
the molecular magnet and the right spin reservoir are 
effectively decoupled in this low-energy subspace. Fig. 
Od)-(f) shows this switching behavior as a function of 
applied electric field. 

In Fig. [2] we assumed in plane Heisenberg exchange in- 
teraction between reservoirs and vertices for all AF reser- 
voirs. This is reasonable, since the AF reservoirs which 
we consider are themself anisotropic. For FM reservoirs 
we assumed isotropic coupling J[ and J' 2 . Additionally, 
we assumed that the third vertex of each molecular mag- 
net is coupled to a separate reservoir by an Ising-like 
interaction with strength Jg = J'. The reason behind 
this assumption is that in this way the sole effect of the 
equilibrium magnetization of the reservoirs on the Hamil- 
tonian of the molecular magnets is to act as an approx- 
imate effective magnetic field J'Se z . This effective field 
for the parameters in Fig. [2] (a)-(g) is on the order of 
75 mT, and can take the role of B', which reduces the 
required external magnetic field. 



V. EXPERIMENTAL REALIZATIONS 



IV. TRANSISTOR BEHAVIOR 

Next, we show how our setup can be used as a logic 
switch whose working is controlled by an external mag- 
netic field. For the first example we assume that ksT <C 
HujABif^B <C Dm-, so that we only need to take the 
states with chirality -1 into account. Referring back 
to Fig. [TJc), we see that spin transport through the 
molecule will be strongly suppressed for magnetic field 
differences such that was Cwb, since in this regime the 
vast majority of the excitations in the reservoirs lack the 
required energy to induce a spin-flip on the molecule; For 
w as i$ wb, transport increases rapidly with was- Hence, 
for magnetic field gradients was ~ w&u our setup can 



A single magnetic dipole moving with constant ve- 
locity v gives rise to a magnetic field Bdi P (r) = 
i^^^pr [3 (e z • e r ) e r — e z ] as well as an electric field 
Edip( r ) = v x B,jip(r) in the laboratory frame. Con- 
ceivably, it is therefore possible to use the setup depicted 
in Fig. G2(a) to measure the switching behavior of a collec- 
tion of spinon-transistors due to the spin-electric coupling 
at temperatures of ~ 1 K. This can be done by measuring 
the difference in voltage drop between points (0, 0, r) and 
(0, r, 0) in the insulating- [at dE' s» Dm/3] and conduct- 
ing (at dE' w 0) state. For the parameters in Fig. [2jh), 
and for AB « 200 mT, the difference in spinon current 
between the two states is ~ 3- 10 10 spinons s^ 1 . For r = 1 
/xm, this leads to a difference in voltage drop of ~ 10~ 13 
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FIG. 3. (a) Proposed setup to measure the switching effect due to the spin-electric effect. The setup consists of a single-layer 
crystal of triangular magnets, of dimension 0.2 x 0.2 /im. The crystal contains ~ 4 • 10 4 molecules, given a lattice constant of 
1 nm. The crystal is weakly exchange-coupled to a bulk collection of AF ID spin chains, such as is realized in SrCu02 (Ref. 
l25l ) or CS2C0CI4 (Ref. [2^. Assuming the lattice constant is commensurate with that of the crystal, the setup contains ~ 4 ■ 10 4 
parallel transistors. The ~ 10 3 molecules at the edge can be accessed electrically by an array of STM tips. Using our previous 
estimate for d and the values from the text, fields between 10 8 — 10 4 V m _1 are required to switch between insulating and 
conducting state, (b) Combining three transistors into a single NAND-gate. (c) Proposed setup to enhance the magnetic field 
due to the accumulated magnons. The crystal of molecular magnets is placed directly above the hole in the superconductor. We 
denote by / the portion of the total flux directly on the magnet that can be enhanced into the area of the molecular magnets. 



V, which is within experimental reachi^i The strength of 
the required switching field E' can be achieved near a 
STM-tip for molecules with reasonable d, see the caption 
of Fig. Ela). This experiment would be interesting in 
its own right, since to our knowledge there have been no 
measurements of the spin-electric effect yet. 

By combing three transistors as shown in Fig. E{b), we 
can create a purely magnetic NAND-gate. The NAND- 
gate is a two-bit gate that gives a logical as outcome if 
and only if both the input bits are 1, and yields a logical 1 
otherwise. We will show how it is possible to implement a 
NAND-gate consisting of magnon-transistors using read- 
ily available materials at ~ 10 K, and we will indicate the 
requirements for a working NAND-gate at room temper- 
ature. A finite non-equilibrium magnetization of a termi- 
nal (reservoir) encodes the logical state 1; the logical state 
has only the equilibrium magnetization present. We 
propose to use the magnetic dipole field B mag = £? mag e 2 
due to the excess N accumulated magnons in the 1-state 
of the terminal as gate-field at the points A-C in Fig. 
OJb) . We use the setup in which each reservoir is coupled 
to all three vertices equally, with strength J'. We assume 
that one reservoir is coupled antiferromagnetically, and 
the other ferromagnetically. In this way, the effects of the 
equilibrium magnetization of the reservoirs on the state 
of the molecular magnet cancel each other. We use the 
fact that for was > w s mag the system is insulating, and 
for was < the system is conducting. Due to ther- 

mal fluctuations, the magnetic field B mag is not constant. 
These fluctuations limit the fidelity of our NAND-gate. 
We characterize the fluctuations by the standard devia- 
tion a = \J ([n — (n)] 2 ) of the number of magnons n on 
the gate, and calculate a using the equilibrium distribu- 



tion of the magnons using the grand canonical ensemble. 
We consider first 2D FM reservoirs with spin S = 10, 
J/ks = 5 K, and lattice spacing a = 1 nm (the ap- 
proximate values for YIG). For the source and drain, we 
consider a single layer sample of dimensions 300 nm x 
300 nm x 1 nm. We put J' /k B = 1 K. As gate, we 
use a crystal of molecular magnets with -Dm/^b = 0.5 
K, and Jm/^b 3> 10 K. Assuming matching lattice con- 
stants, this setup contains 300 parallel single-molecule 
transistors. We estimate the field at the position of the 
molecular magnet [see Fig. E^c)] due to the n magnons as 
|B mag | « \^ng hb f / {da x a y ) . To get rid of the magnetic 
field due to the equilibrium magnetization, we can add 
an additional FM with antiparallel equilibrium magneti- 
zation and no excess magnons on the opposite side of the 
molecule. Additionally, we use parameters gu = 9r = 2, 
/ = 0.15, B = 200 mT, B' = 60 mT, and AB = 50 mT. 
In the conducting state, we find a magnon current ex- 
ceeding 1.5- 10 10 s _1 , which amounts to a switching time 
of a transistor of ~ 300 ns. In our model, the fidelity of 
a single-NAND-gate then exceeds 99.9%. 

By using materials with an increased g- factor, we can 
create a NAND-gate that functions at room temperature, 
consisting of transistors with a ~ 11 ns switching time. 
We use the same setup as in the previous paragraph, but 
with parameters S = 3, J/ks — 500 K, a — 1 nm for the 
2D FM reservoir; D M /k B = 2.5 K, and J M /k B > 300 
K for the molecular magnet; and J'/fcs = 100 K, <?m = 
gn = 20, / = 0.15, B Q = 1 T (easily achievable near 
the surface of a FM), B' = 100 mT, and AB = 50 mT. 
We find a magnon current exceeding 2.7 • 10 11 s _1 in the 
conducting state, which corresponds to a 11 ns switching 
time. Reducing the magnon fluctuations on the gate can 



7 



further reduce the switching time. As before, the fidelity 
of a single NAND-gate exceeds 99.9%. 



VI. DISCUSSIONS 

In this section, we will discuss several different require- 
ments that have to be fulfilled for our perturbative cal- 
culations of the tunneling current through the molecular 
magnet to be valid. The first constraint on our calcula- 
tions concerns the validity of our spin wave analysis of 
the FM spin chains; the number of magnons per site has 
to satisfy (afa^) <C 2S. In our calculations, the aver- 
age number of magnons per site is typically 0.05-0.13, 
so that non-interaction spin wave theory is valid. The 
AF theory is valid for energies much smaller than the 
exchange interaction J. 

We checked the validity of our sequential tunneling ap- 
proach in a self-consistent manner. For the AF reservoirs, 
the criterion is simply that the tunneling current is much 
smaller than the current in the ballistic system, that is 
Is(AB) <C guB^B jh. For the FM reservoirs, we require 
that the broadening of the energy levels of the molecule 
is smaller than the unperturbed level splitting. In other 
words, all transition rates R^J R of the FM system satisfy 

^^/riB^-^^^-ejl/h, (13) 

where e°, e° are the unperturbed energies of the states 
l*>> I/)- 

Another constraint is given by the fact that, near Breit- 
Wigner resonances, the current through the molecular 
magnet can be strongly increased due to coherent tun- 
neling processes. This only holds at low temperatures, at 
higher temperatures the broadening of the thermal dis- 
tribution destroys coherent tunneling, and the sequential 



tunneling approach is valid again. The minimal temper- 
ature T for FM reservoirs is given by 

Rf/ R /n B (LJi - u f ) < k B T/h. (14) 

All our calculations are at high enough temperature for 
the sequential tunneling approach to be valid for the FM 
system. 

We note that relaxation of the state of the molec- 
ular magnet can be neglected as long as the coupling 
strength between the reservoirs and the molecular mag- 
net exceeds the coupling between the molecular magnet 
and hyperfine- and phonon baths. 

VII. CONCLUSIONS 

Using a sequential tunneling approach, we have studied 
transport of magnons and spinons through a triangular 
molecular magnet which is weakly coupled to two spin 
reservoirs. We have shown that, by changing the state of 
the molecular magnet through application of an electric- 
or magnetic field, we can control the magnitude of the 
spin current through the molecular magnet. We used 
this fact to propose a magnon-transistor, whose operation 
can be controlled by an electric- or magnetic field. We 
have shown for which parameters our transistor could 
operate at room temperature with a 11 ns switching 
time. We have shown how several magnon-transistors 
can be combined to create a NAND-gate. 
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